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Lambda production in the DIS target fragmentation region 
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IFPA, Universite de Liege, B4OOO, Liege, Belgium 

By using a recently obtained set of Lambda fracture functions, we present predictions 
for Lambda production in the target fragmentation region of Semi-Inclusive Deep 
Inelastic Scattering in CLAS@12 GeV kinematics, supplemented with a conservative 
error estimates. We discuss a number of observables sensitive to the assumptions of 
the underlying theory and many of the assumptions of the proposed phenomenolog¬ 
ical model. 
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I. INTRODUCTION 

Hadron production in Semi-Inclusive Deep Inelastic Scattering (SIDIS) is usually de¬ 
scribed in terms of universal parton distributions and fragmentation functions. Thanks to 
the factorisation theorem, hadronic cross sections are obtained by convoluting short-distance 
partonic cross sections, calculable in pertubation theory, with such distributions. To lowest 
order in the strong coupling, this mechanism is expected to describe hadron production 
in the so called current fragmentation region, i.e. the phase space region in which the 
struck parton hadronises. In order to obtain a global description of the particle production 
spectrum, and in particular of hadron production in the target fragmentation region, the in¬ 
troduction of new non-perturbative distributions is mandatory. This issue was early realised 
in Ref. [l| where the concept of fracture functions was introduced. The latter parametrise 
the hadronisation into the hnal-state hadron of the coloured spectator system which results 
from the removal of the scattered parton from the initial-state hadron. For this reason, their 
flavour and energy dependecies are expected to be significantly different from fragmentation 
functions which parametrise the fragmentation of a single parton into the observed hadron. 
Fracture functions, by construction, simultaneously encode information both on the parton 
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partecipating the hard scattering and on the fragmentation of the spectator system into the 
observed hadron. Therefore they constitnte the connection between forward particle pro¬ 
duction at small transverse momentum (i.e. target fragments) and high momentum transfer 
processes (i.e. DIS). Although intrinsically of non-perturbative nature, the scale dependence 


of such distributions can be predicted by perturbative QCD [l|. Fracture functions obey, 


in fact, DGLAP inhomogeneous evolution equations which result from the structure of 
collinear singularities in the target fragmentation region [l, Moreover, a dedicated fac¬ 


torisation theorem j^, 5| guarantees that fracture functions are universal distributions, at 
least in the context of SIDIS. Among baryons. Lambda hyperons are predominantly pro¬ 
duced in the SIDIS target fragmentation region and show a signihcant leading particle effect, 
i.e. they carry a signihcant fraction of the incoming proton momentum. For such reasons 
they have been used as a case study and a hrst attempt to determine Lambda fracture 
functions have been recently presented in Ref. Q| by performing a simultaneous QCD £t to 
a variety of Semi-Inclusive Lambda production data collected in lepton-nucleon scattering. 
In the present paper, by using this model, we present predictions for Lambda observables 
in the target fragmentation region of neutral current (NC) Deep Inelastic Scattering (DIS) 
focusing on CLAS@12 GeV kinematics. The paper is organised as follows. In Sec. [TTl and 
Sec. uni we hrst briehy review the Inclusive and Semi-Inclusive DIS cross sections in lepton- 
nucleon scattering. In Sec. [IV]we review some details of the modelisation of Lambda Fracture 
Functions. In Sec. El we present and discuss a number of observables sensitive to the as¬ 
sumptions adopted in the model which can be used to further constrain it. Finally in Sec. I VI I 
we summarise our results. 


II. INCLUSIVE DIS 


The deep inelastic scattering cross section of a lepton I oh a proton p 

l{k) + p{P) ^ l{k') + X , 


( 1 ) 


with four-momenta k and P, respectively, is usually described in terms of the invariants: 




P ■ q 




_ V _ ^ ^ /n2 _ 2 /o^ 

'"^~2P-q' P-k~ {s-ml)x' 

where k' and q = k — k' are the outgoing lepton and virtual boson four-momenta, respectively, 
s = {P + k)‘^ is the centre of mass energy squared and = sy{l — x) + rri^ is the invariant 
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FIG. 1: Kinematic coverage at CLAS@12 GeV after cuts in eq. m- 


mass squared of the hadronic final state, with rrip the proton mass. The leading order NC 
DIS cross section for the scattering of an electron of energy Ee on a proton target then reads 

(f(jep^eX 2Tial 


dxsdQ'^ 




■(1 + (1 - yf) XI [fq/pi^B, Q^) + fq/p{XB, Q^)] , 


(3) 


where the sum runs over active quarks q with electric charge Cq. The differential cross section 
in eq. ® is evaluated by using free-nucleon, leading order, parton distributions fq/p{xB, 
of Ref. [7|, setting the factorisation scale to y\ = . In order to provide a minimal quark- 

flavour separation, we consider both electron scattering on proton and deuteron targets. 
Cross sections on the latter are obtained by averaging cross sections on proton and neutron 
targets. The latter are obtained applying isospin symmetry, i.e. by exchanging u -H- d and 
it -H- d parton distributions in eq. ([3]). We set the electron beam energy Eg to 12 GeV. The 
label G stands for a set of temptative cuts which dehne the NC DIS selection: 


0.2 <y< 0.8, 1 GeV^ > 5 GeV^. 


(4) 


The resulting phase space coverage is shown in the {xbtQ"^) plane in Fig. ([T]). In order to 
reduce the dependences on higher order corrections, all predictions presented in the following 
are normalised, if not otherwise stated, to the inclusive NC DIS cross section, which is 
obtained integrating eq. ([3]) over the phase space region D dehned by constraints in eq. (jlj). 
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III. SEMI-INCLUSIVE DIS 


We consider the Semi-Inclusive process 


l{k) + p{P) ^ l{k') + k{h) + X , 


( 5 ) 


where, beside the scattered lepton, an additional Lambda hyperon is detected in the hnal 
state with four-momentum h. Final-state hadrons produced in SIDIS are generally described 
by using the Lorentz-invariant variable 


P-h 


Et 


Zh = 


1 -I- cos 9* 


( 6 ) 


P-q E’p(1-xb) 2 

The last equality holds in the photon-hadron centre-of-mass frame, with the photon mo¬ 
mentum aligned in the +z direction and 9* the hadron production angle with respect to 
the photon direction. Hadrons produced collinearly to the spectator system have 9* ~ vr, so 
that, in terms of the Zh variable dehned in eq. (E]), they overlap with soft ones (for which 
instead ~ 0 irrespective of the production angle) and both accumulate at Zh — 0. The 
Zh variable dehned in eq. (E]) is therefore well suited to described hadron production in the 
current region, but presents an ambiguity in dealing with hadrons produced by target frag¬ 
mentation. In order to avoid this problem, cross sections can be evaluated in terms of the 
energy fraction zg dehned by 

E 


Zg = 


K 


2E* 


E*p{1-xb) IF’ 


c = — 


(7) 


where E*{1 — xb) = W/2 is the spectator energy in the photon-hadron centre-of-mass frame. 
Adopting such dehnition, higher order corrections can be systematically taken into account, 
both in the current and in the target fragmentation region so that diherent hadron 
production mechanisms are distinguished only by their peculiar zg spectrum. Adopting 
these dehnitions, the neutral-current semi-inclusive lowest order cross section for producing 
an unpolarised Lambda oh a proton in the target fragmentation region reads 

fj^ ^-- 

= -^-^(1 + (1 - + M,i,(xb,C,,Q")] . ( 8 ) 

The cross section has been expressed for later convenience in term of the C, variable 
{xb + C < 1) ill eq. ([7]), and the Jacobian J = (■[(! — has been explicitely in¬ 

dicated j^. The latter reduces to unity in the high-energy limit and it is therefore often 








5 


omitted in the literature. In eq. dS]) the production of unpolarised Lambdas in the remnant 


Q. 


These distributions express 


direction is described by fracture functions 
the probability to hnd a parton of flavour i with fractional momentum xs at virtuality /i|^ in 
the proton conditional to the detection of a target Lambda with a fraction C, of the incoming 
proton momentum. As for inclusive parton distributions, we set the factorisation scale to 
hF = order to obtain cross sections on isoscalar target we proceed as in the inclusive 

DIS case, exploiting isospin symmetry of the initial conditions. More details on this point 
may found in Sec. lIVi Distributions in a given kinematic variable v = v{xb, Q‘^,C) ci-re then 
calculated integrating the SIDIS cross section in eq. (j8]) as follows 

1 


Avi Avj 


'’^dxsdQ^<K 


©(■y - Vi) 9(vi+i - v ), 


(9) 


where the index i labels the Lth bin, Uj+i and Vi indicate the experimental bin-edges, Avi = 
Vi+i — Vi stands for the bin-size and D' is a subset of the DIS selection, D' C D. 


IV. MODEL DETAILS 


Lambda fracture functions appearing in eq. ([8]) have been determined through a global 
QCD £t to a variety of semi-inclusive Lambda production data in Ref. jb]. In that analysis 
we assumed that, at an arbitrary low but still perturbative scale Ql, fracture functions 
factorise into the product of ordinary parton distributions fi/p{xB-, Ql) and what we address 
as spectator-fragmentation functions D^^^{zg)- 

(1 - Xb) {xb, C, Ql) = M^ipixB, Zg, Ql) = fi/p{xB, QI)D^/p{zg) , i = q,q,g. (10) 

Such an assumption, supported by the £t, is motivated by considering the relevant timescales 
in the process. The hard scattering, controlled by parton distributions, occurs in fact on 
timescales 0{1/Qq) much shorter than the typical timescale of the fragmentation process, 
0{1/ Aqgd), controlled by spectator fragmentation functions. These initial conditions for 
fracture functions at are then evolved to scales relevant for the experiments and the 
parameters controlling D^^^{zg) extracted by performing a £t to available data. The scale 
Ql is fixed in the fit to 0.5 GeV^. As discussed in Ref. |6|, the latter were essentially able 
to constrain only a subset of the initial conditions in eq. ffTOj) . in particular the spectator 
fragmentation functions of valence u and d quarks and that of sea-quarks, for which a. 
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FIG. 2: Spectator fragmentation functions D^^^{zg) for valence u (green), d (red) and sea (yellow) 
quarks extracted from the fit of Ref. IqJ. The bands represent the propagation of experimental 


uncertainties according to the = 1 criterion. 


flavour-independent, common function was assumed: 

Mt/p{xB,ZG,Ql) = fu^/p{xB,Ql)D^^/p{zG) , 

Mi/p{xB,ZG,Ql) = fd,/p{xB,Ql)Dl/pizG) , ( 11 ) 

ZGiQq) = fqB/pi^BiQo)Dq^^p[zG), Qs — Ug, dg^Ug, . 

Such distributions are shown in Fig. ([2]) as a function of zg'- the fragmenting spectrum 
of the ud-spectator into Lambdas, is harder with respect to the uu-one, The 

leftover spectator system, associated with the hard scattering occurring on sea quarks, 

has higher Fock components and shows a softer spectrum. The £t was insensitive to any 
reasonable choice of gluon spectator fragmentation function, so the latter was hxed to be 
equal to the sea-quark one, D^Jp^zg) = D^^Jp^zg). The structure of the initial conditions 
allows to evaluate valence-quark fracture functions at any as the difference between 
M^Ip and The cross sections on a neutron target requires the knowledge of neutron- 

to-Lambda fracture functions, We relate the latter to assuming the following 
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relations: 


Miln{xB:ZG,Q‘^) 

Mq,jn{xB,ZG,Q‘^) 


Mtlp{.XB,ZG,Q^), 

Ml/pixB,ZG,Q'^) , 

Mg^/pi^B^ZG^Q"^), q's=Us,ds,... qs = ds,Us, 


( 12 ) 


The first one in eq. flT^ appears to be natural as the quark content (ud) of the leftover 
spectator system is the same on the left and right hand side. The second assumption 
implies that the (id-spectator has the same fragmenting spectrum as the uu-one and it 
does not have any physical motivation other then reducing the number of free parameters 
in the £t. The third one relies on the particular flavour-symmetric choice for DgA^ZG) 


in eq. m- Following the method outlined in Refs, j^, [l^ we provided, beside the best 


fit parametrisation, additional 14 Lambda fracture functions alternative parametrisations 
satisfying the = 1 criterion. In this way experimental uncertainties can be propagated 
to any other observable by computing it for each given alternative set and then adding in 
quadrature the displacements with respect to best fit result. This method has been used 
to obtain the error bands associated with the spectator fragmentation functions in Fig. ([2]) 
and will be used in the following to estimate experimental uncertainties for the relevant 
cross sections. We close this Section mentioning that, at least in principle, we would be 
interested in promptly produced Lambdas. It is well known, however, that a fraction of 
the measured Lambda yield comes from the decay of heavier resonance into Lambdas, the 
so-called feed-down effect. The subtraction of these fractions from the Lambda yields was 
not clearly stated or even technically achievable in many of the experimental analyses whose 
data have been used in the fit of Ref. jb]. We assumed therefore that the quoted yields 
referred to an unsubtracted Lambda sample. These production mechanisms are taken into 
account by the Lambda fracture function set of Ref. jb], via effective modifications of the 
spectator fragmentation functions returned by the fit. Such an assumption must be kept in 
mind when comparing predictions based on the present model with forthcoming data. 


V. PREDICTIONS 

Cross section differential in the energy ratio zg defined in eq. ([7]) characterises the full 
particle production spectrum. The latter is given by the sum of the target fragmentation 






contribution, given in eq. ([H]), and the current one, in which fracture functions appearing 
in eq. (|H]) are replaced by appropriate products of parton distribution and fragmentation 
functions. As discussed in Section m different hadron production mechanisms are then 
distinguished only by their peculiar zg spectrum and target fragmentation can be quantihed 
without imposing any arbitrary kinematical cuts. So far, unfortunately, experimental data 
have not been presented in terms of this variable. Even in that case, however, the extraction 
of fracture functions with such a procedure requires an accurate knowledge of the current 
fragmentation contribution at low scales, whereas fragmentation functions are generally 
constrained at much higher scales than the ones involved in SIDIS experiments. In order 
to circumvent this problem we assumed in Ref. 6j that current and target fragmentation 
give their dominant contributions in distinct regions of space phase. Within this context it 
proves useful to introduce the Feynman’s variable 


Xp = ±{z, 


'G 


- 


(13) 


dehned in the photon-proton centre-of-mass frame. The parameter e will be used in the 
following to estimate the sensitivity of the predictions to Lambda mass corrections and it is 
hxed to e = 1. We assumed that the current and target contributions can be kinematically 
separated in terms of this variable with target fragmentation giving its contribution for 
xp < 0 and current fragmentation for > 0. We stress again that the choice of = 0 
as a sharp separation point is arbitrary and frame dependent. Moreover it is reasonable to 
expect that there will be an overlap region in which both fragmentation mechanisms will 
contribute, as suggested by the 0{as) calculation of Ref. j^. This strategy was adopted 
in Ref. jb] in the extraction of Lambda fracture functions and we shall consider it as an 
operative choice to be tested against forthcoming data. We have introduced in eq. flT3|l the 
Lambda transverse mass, rrip = Pa ± + dehned in terms of its transverse momentum 
and mass squared. Since it is experimentally known ll| that {p\±) -C m\, we approximate 
~ m\ with tua = 1115.683 MeV 12|. Lambda-mass effects, introduced via eq. flT5]l . are 
sizeable at low energies and moreover not compatible with the pQCD factorisation theorem. 
As described in Ref. [l3|, such corrections are applied to the Lambda leptoproduction cross 
sections in eq. (E]) via the extra J factor. We begin our overview of results presenting 
in Tab. ([T]) predicted yields and absolute cross sections for the production of Lambdas with 
Xi? < 0 within the DIS selection dehned in eq. dH). Results on proton and deuteron targets 
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Target / Observable 

(n(A)) 

proton 

deuteron 

0.038 ± 0.003(exp)lg QQ 4 (mass)^QQQ^(scaA) 
0.032 ± 0.002(exp)lgQQ4(mass)^QQQj;(sca/e) 

Target / Observable 

[pb] 

proton 

deuteron 

2382 ± 17Q{exp)f(^g{mass)f)^^{scale) 
1758 ± lt)2{exp)f}§^{mass)f]^{scale) 


TABLE I: Predicted Lambda yields and cross section for xp < 0 on proton and deuteron targets. 


Quoted errors represent the propagation of experimental uncertanties from the /li (exp), the 
expexted sensitivities to mass corrections (mass) and to factorisation scale variations (scale). 


are shown. The former are larger than the latter, since, as already mentioned, onr model 
returns larger and harder than Both the yields and the absolute cross sections 

are supplemented by errors. The hrst one corresponds to the propagation of experimental 
uncertanties coming from the £t jb] and it is denoted with the label (exp). It amounts 
to an average uncertainty of 7% on the yields. Among theoretical errors we address the 
sensitivity to Lambda mass corrections and higher order corrections. The former is assessed 
by varying, arbitrarily, the parameter e appearing in eq. flT3|) in the range e G [0.9,1.1] and it 
is indicated with the label (mass) in Tab. ([T]). The latter is assessed, as costumary, by varying 
the factorisation scale both in fracture and ordinary parton distributions in the range 
p.% G [0.5 Q^, 2 and it is indicated with the label (scale) in Tab. (jl]). From these numbers it 
appears that there is, given the relatively low beam energy of the experiment, a rather large 
sensitivity to Lambda mass corrections. On the other hand, errors associated to estimated 
higher order corrections are smaller than experimental uncertainties. The yields appears to 
be particularly stable against scale variations since the factorisation scale is simultaneously 
varied both in the numerator (i.e. fracture functions) and denominator ( i.e. inclusive 
parton distributions). We stress here that while the quoted experimental uncertainties have 
a precise statistical meaning the other two must be considered as temptative estimations 
of systematics errors associated with theoretical predicitions. We present in Fig. ([3]) the 
normalised xp spectrum supplemented by the correponding uncertanties. We now turn 
our attention to the production properties of backward Lambdas. The averaged values 
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FIG. 3: Normalised single-differential cross sections as a function of xp- Cross sections on pro¬ 
ton (H) and deuteron (D) targets are shown. From top to bottom the error bands represent the 


propagation of experimental uncertainties from the fit 
factorisation scale variations. 


a/, sensitivity to mass corrections and to 


of kinematical variables for Lambda production in NC DIS with < 0 are given by 
{xb) = 0.17 (0.18), fy) = 0.48 (0.45), (Q^) = 2.0 (1.9) GeV^ and = n.i ( 10 . 3 ) GeVl 
In parenthesis we have indicated the corresponding values for the inclusive DIS case. The 
values for the two processes are quite close to each other, a feature which should be ascribed 
to the factorised ansatz for fracture functions at the lowest scale in eq. flT 0 |l which is almost 
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xb y 



FIG. 4: Normalised single-differential eross seetions as a function of leptonie kinematic variables 
with the additional requirement of deteeting a Lambda with xp <0. Cross sections on proton and 
deuteron targets are shown. 

preserved by evolution from Qq to (Q^). In Fig. (jl]) we present the normalised single¬ 
differential cross sections as a function of scattered lepton variables with the additional 
requirement of detecting a Lambda in target region, xp < 0. Such distributions show a 
qualitatively similar shape irrespective of the target considered, either proton or deuteron. 
The Xb distributions peaks at the lowest accessible values of xb, given by the boundary in 
Fig. ([H). The bulk of the cross section resedes, as expected, at very low values of and 
invariants. The combined study of such distributions together with the corresponding 
ones in inclusive DIS can potentially highlight correlations between the hard scattering and 
the spectator fragmentation into target Lambdas. In Fig. (jS]) we present the Lambda single¬ 
differential cross section as a function of xp, integrated in ranges of W^. The very backward 
production regime (at large and negative xp) is accessed only at highest values of W^. This 
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FIG. 5: Normalised Lambda single-differential cross section as a function of xp integrated in 
different range ofW‘^ on proton (H) and deuteron (D) targets. 


is the combined effect of hadron mass corrections, via eq. ca. and the energy spectrum 
of the spectator fragmentation functions shown in Fig. ([2]). In the top left panel of 
Fig. ([6]) we show the xp dependence of the cross section in ranges of Q^. As already seen 
in Fig. dl]), the bulk of the cross section is at low Q^, althought it remains non-negiglible to 
the highest accessible Q^. The measurement of the Semi-Inclusive cross section in the target 
region as a function of can validate the leading twist nature of particle production in this 
region of phase space, as assumed by fracture functions formalism. In the top right panel 
of the same hgure we present Lambda multiplicities in ranges of as a function of xp. In 
this case distributions are normalised to where the additional index i = 1,2,3 stands 
for the corresponding range indicated on the plot which supplements the DIS selection 
n. A mild rise of the multiplicity can be observed as increases, which can be possibly 
ascribed to the QCD evolution of fracture functions. It would be extremely interesting to 
compare these distributions with the corresponding one in photoproduction regime in order 
to determine to which extent the transition to the non-perturbative regime in affects the 
Lambda spectrum in the target region. In the bottom row of Fig. ([H]) we show the quark- 
flavour decomposition of the single differential cross sections as a function of xp normalised 
to cTqJ®, where hi' stands for the DIS selection D supplemented with the cut xp > 0.25. 
Assuming that the experiment can be performed both on proton and deuteron targets, the 
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FIG. 6: Top left: Normalised Lambda single differential cross section as a function of xp in different 
range of on a proton target. Top right: Lambda multiplicities as a function of xp in ranges 
of Q^. Bottom: quark-flavour decomposition of the Lambda single-differential cross section as a 
function of xp on a proton (left) and deuteron (right) target with the additional cut xp > 0.25 
imposed. 

plots show that such a selection can provide an optimal valence quark-flavour discrimination 
for Lambda fracture functions. In the proton target case the Lambda spectrum is dominated 
by scattering on valence n-quarks with maximal sensitivity to the fragmentation of the ud- 
spectator system into Lambdas. In the deuteron case, the generalised isospin relations in 
eqs. flT^ allow the extraction of the uu-spectator fragmentation functions. In Fig. ([7]) the 
normalised single differential cross sections as a function of xp are presented in and 
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1 < 0^ < 2 GeV^ 2 < 0^ < 4 GeV^ 4 < < 16 GeV^ 



FIG. 7; Normalised single differential cross sections as a function ofxp, l/aQ^da^/dxp, in various 
bins of xb and Q^. Distributions for hydrogen and deuteron targets are presented and correspond¬ 
ingly normalised. 

bins. This way of presenting the data is probably the more exhaustive and it might 
be valuable for the determination of Lambda fracture functions in forthcoming global fit 
analyses. We conclude this Section presenting in Fig. ([8]) the normalised Lambda single¬ 
differential cross sections on a proton target as a function of in four different bins of 
Xb- The final-state Lambda is required to have Xii’ < 0. The Q^-differential cross section 
deserves special attention since this observable may provide crucial test for the predicted 
evolution of fracture functions and validate the key assumptions of the underlying theory. 
Given the relatively low values of accessed by the experiment, the spectrum shows 
significant hadron mass corrections, as can be inferred comparing default predictions with 
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Q" [GeV^] [GeV^] 


FIG. 8: Normalised Lambda single-differential cross section on a proton target as a function of 
in different bins of xb and integrated in the range xp <0. 

the one in which the Lambda mass has been set to zero. Their effect is to suppress the cross 
section as xb increases. In view of these results, the genuine dependence of the cross 
section and mild logarithmic effects generated by QCD evolution of fracture functions can 
get obscured by hadron mass corrections. Therefore the interpretation of forthcoming data 
will require a proper modelisation of the latter either with the basic method described in 
this paper or with more refined treatment as the one discussed in Ref. jl^ . 


VI. CONCLUSIONS 

In this paper we have considered Lambda production in the target fragmentation region 
of electron-proton deep inelastic scattering. We have presented, based on a recently ob¬ 
tained set of Lambda fracture functions, predictions for a number of relevant observables 
supplemented with a conservative error estimates. In future perspective, the subdivision of 
xp spectra in bins of and xp can be a valuable input for forthcoming hts. Given the 
energy range of the considered experiment, the possibility to use different light targets offer 
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an additional handle on Lambda fracture functions quark-flavour separation in the valence 
region. The study of the of the dependence of the cross sections can be valuable to test 
and validate the key feature of the underlying theory and many of the assumptions of the 
proposed phenomenological model. 
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